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1.0  INTRODUCTION 


A  thorough  understanding  of  target  backgrounds  is  necessary  to  describe  the  perform¬ 
ance  of  passive  infrared  (fRf  surveillance  systems.  A  model  has  been  developed  based  on  the 
well  known  Cox-Munk  wave-slope  statistics  (Cox «nd  Munk,  1954)  and  the  LOWTRAN  6 
optical  propagation  code  (Kneizys  et  al.,  1983),  to  predict  the  behavipUr  of  sea  radiance  back¬ 
grounds.  Similar  approaches  have  been  published  in  the  past,  but  they  neglected  important 
processes.  Saunders  (1968)  modeled  sea  radiance  but  neglected  the  path  radiance  between  the 
sensor  and  the  sea  surface  and  also  aerosol  effects.  Sidran  (1981),  as  well  as  Schwartz  and  Hon 
(1986)  studied  only  the  effects  of  sea-surface  roughness  on  emissivity.  Wilf  and  Manor  (1984) 
addressed  the  path  radiance  between  the  sensor  and  the  sea  in  a  simplified  way  by  neglecting 
aerosols.  The  model  presented  here  is  an  improvement  over  earlier  attempts  since  it  addresses 
the  effects  of  refraction  and  molecular  absorption  and  emission,  as  well  as  the  aerosol  effects, 
which  are  important  in  determining  the  path  radiance.  The  reflection  geometry  and  the  wave- 
slope  distribution  determine  those  portions  of  the  sky  that  reflect  energy  into  the  line-of-sight 
of  the  receiving  system.  A  provision  is  included  for  determining  the  effects  of  wave  shadowing. 
The  surface  reflectivity  and  emissivity  are  calculated  using  the  complex  index  of  refraction  of 
sea  water  in  Fresnel’s  formulas.  The  model’s  ability  to  describe  the  sea  and  sky  radiances  for 
zenith  angles  close  to  the  horizon  is  demonstrated  by  comparing  calculated  radiances  using 
radiosonde  soundings  of  the  atmosphere  with  measurements  of  the  8-  to  12-^m  radiances  using 
a  calibrated  thermal-imaging  system. 


2,0  THEORY 


2.1  BASIC  IDEAS 

A  IR  system  looking  down  at  the  .sea  surface  receives  energy  from  different  sources. 

The  air  between  sensor  and  sea  surface  emits  and  scatters  radiation  into  the  line-of-sight  of  the 
system.  The  sea  surface  reflects  parts  of  the  sky’s  radiance  and  emits  with  the  water  tempera¬ 
ture  T^.  The  emissivity  and  reflectivity  of  salt  water  can  be  determined  by  Fresnel’s  formulas. 

A  frequency-dependent  complex  index  of  refraction  has  to  be  used.  Reflected  sky  radiance,  as 
well  as  the  thermal  emission  of  the  sea,  have  to  be  multiplied  by  the  transmissivity  t.  of  the 
path  between  sensor  and  sea  surface  to  address  loss  of  energy  due  to  scattering  and  absorption. 
Figure  I  shows  the  geometry  of  a  reflection  due  to  waves  with  an  angular  tilt  of  a  in  wind 
direction  and  P  perpendicular  to  the  wind  direction.  The  curvature  of  the  path  was  a  result  of 
atmospheric  refraction. 

Equation  (1)  determines  the  radiance  received  by  the  sensor  and  contains  all  the  above 
mentioned  contributions. 

’’a^o) 

S(0,k)  =  J"  r^frg)  (1  TrlBBfk.Tglr^))  dr^  (1) 

T-afS,) 

+  7(s,)  R(n,k)  J  (|,,,r)BB(k.Ta(7a))  dr^ 

7-a(S,) 

+  T(S|)  (1  -  R(ak))  (l/7r)BB{k,T^)  . 
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Figure  1.  The  geometry  of  a  reflection  at  a  wave  is  shown.  The  sensor  zenith  angle  0  is  modified 
by  atmospheric  refraction  resulting  in  an  zenith  angle  0.  The  reflection  angle  and  the  zenith 
angle  ^  from  where  radiance  is  reflected  in  the  sensor  line-of-sight  are  defined  by  the  wave  slopes 
a,  (3,  and  the  angle  0. 


To  evaluate  equation  (1),  determine  the  transmission  the  optical  path  s.  and  the  reflec¬ 
tivity  R.  The  LOWTRAN  6  model  is  used  to  determine  these  values  numerically.  The  actual 
optical  path  is  controlled  by  the  wave  tilts  Sx  and  Sy  (Sx  =  tan  a,  Sy  =  tan  )8),  the  sensor  zenith 
angle  0,  and  refractivity  effects  (i.e.,  the  vertical  distribution  of  temperature  T^,  pressure  (p), 
and  relative  humidity  (RH)).  Instead  of  momentary  values  of  radiance,  we  are  looking  at 
spacial  averages.  The  statistical  distribution  of  wave  tilts,  together  with  the  reflection  geome¬ 
try,  determine  the  statistical  distribution  of  zenith  angles  from  radiance  reflected  into  the 
observer  line-of-sight.  The  average  radiance  is  calculated  by  integrating  over  all  zenith  angles 
(equation  2)). 

ra(Si) 

S(0,k)  =  J"  rj(rg)  (1 'rrlBBfk.Tgfrg))  dr^  (2) 


The  wave-slope  distribution  necessary  for  these  calculations  has  to  be  obtained  from  a 
sea-state  model. 
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2.2  SEA-STATE  MODEL 


The  Cox-Munk  statistics  were  used  to  describe  the  wave  tilts.  The  probability  density 
function  of  wave  tilts  for  a  wind-driven  sea  was  given  by  a  Gram  Charlier  series.  The  coeffi¬ 
cients  were  determined  experimentally  by  Cox  and  Munk.  A  two-dimensional  gaussian  distri¬ 
bution  has  been  used  in  this  paper.  Equation  (3)  shows  the  gaussian  function  together  with  the 
equations  to  determine  the  variances  (x  pointing  downwind,  y  crosswind)  of  the  wave 

slopes  as  functions  of  the  windspeed  (WS)  (in  m/s). 

P(Sx,  Sy)  =  l/(27raxOy)  EXP (O-SCS^^/ 0^,2  +  Sy^/oy^)) 

-  0.003  +  1.92  10'^  WS 

ay2  =  3.16  10  ^  WS  .  (3) 

One  perspective  of  this  investigation  was  to  study  how  well  a  simple  sea-state  model 
corresponds  with  measured  radiances.  Different  authors  (Schwartz  and  Hon.  1986;  Wu,  1972) 
showed  that  the  Cox-Munk  model  underestimated  the  total  variance  of  the  wave  slopes. 

Future  work  should  investigate  the  influence  of  different  relationships  between  the  variances 
of  wave  slopes  and  windspeed.  or  other  versions  of  a  sea-state  model. 

2.3  GEOMETRY 

Two  different  coordinate  systems  were  used.  The  wind  oriented  system  with  its  x-axis 
pointing  downwind,  the  z-axis  vertically  upward,  and  the  y-axis  perpendicular  to  x  and  z  so 
that  a  right-handed  system  was  formed.  In  the  wave-oriented  system  the  x'-axis  was  within  the 
wave  surface  pointing  downwind.  The  z'-axis  is  normal  to  the  wave  surface  and  the  y'-axis  is 

perpendicular  to  x'  and  z'  so  that,  again,  a  righthanded  system  was  formed.  The  transition  from 

one  system  to  the  other  was  possible  through  the  use  of  the  wave  slopes  S^.  Sy  given  in  the 
wind-oriented  system.  By  calculating  the  coordinate  transforms  between  both  systems,  it  was 
possible  to  derive  the  zenith  angle  n  (for  those  portions  of  the  sky  from  where  radiance  is 
reflected  into  the  line-of-sight)  and  the  reflection  angle  Cl  (equation  (4)). 

cos  n  -  2S^  A  cos  0'  cos  4>  +  2Sy  A  cos  0'  cos  4’  B  A  sin  0’ 

cos  fl  =  S„  A  cos  0'  cos  4>  +  S,.  A  cos  0'  cos  4>  +  1  A  sin  0' 

X  y 

A  =  5^2  +  Sy2  +  1 

B  =  5^2  +  Sy2  +  1  .  (4) 

2.4  SHADOWING  OF  WAVES 

For  a  given  configuration  of  sensor  zenith  angle  and  wave  slope,  a  certain  percentage 
Pq  of  waves  were  hidden  behind  other  waves.  Thus,  for  this  percentage  of  cases  the  reflected 
sky  radiance  did  not  come  from  the  zenith  angle  that  results  from  the  geometry,  but  from  other 
zenith  angles.  In  contrast  to  another  paper  (Saunders,  1968).  it  is  assumed  that  the  radiance 
reflected  from  those  zenith  angles  also  contributed  to  the  received  radiance  of  the  sensor. 
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In  cases  where  shadowing  must  occur  (0  >  90°  and  fl  >  90°),  it  is  assumed  that  the 
radiance  was  reflected  by  the  other  side  of  the  same  wave.  Equation  (5)  states  zenith  and  reflec¬ 
tion  angle  for  those  cases. 

cos  n  =  -2Sjj/A  cos  0'  cos  $  -  2Sy/A  cos  0'  cos  <I»  -  B/A  sin  0'  and  (5a) 
cos  ft  =  -Sjj/  A  cos  0'  cos  ^  -  Syl  A  cos  0'  cos  4>  +  1  /  A  sin  0'  (5b) 

A  =  +  Sy2  +  1 

B  =  5^2  +  s^2  _  , 

In  cases  of  multiple  reflection  as  well  as  other  shadowing  events,  it  is  assumed  that  the 
shadowed  optical  path  hits  another  part  of  a  normally  distributed  sea  surface.  This  situation  is 
described  by  the  zenith-angle  distribution  already  calculated.  The  probability  for  shadowing  is 
distributed  over  the  zenith  angles  according  to  the  distribution  function  (Appendix  I). 

2.5  REFLECTIVITY  AND  EMISSIVITY  OF  SEA  WATER 

The  reflectivity  (and  emissivity)  of  water  as  a  function  of  the  reflection  angle  ft  and  the 
complex  index  of  refraction  (mr,  mi)  is  given  by  Fresnel’s  equations  (6),  (Stratton,  1941). 

1)2  =  0.5(m^2  m.2  _  +  ((m^^  j^,2  sin2ft)2  +  4m^2rf,.2^0.5) 

v2  =  0.5(-<m^2  _  nij2  -  sin2ft)  +  ((m^2  _  ^,.2  _  sin2fj)2  +  4j^^2^,2)0.5j 

C  =  (cos  ft  -  U)2  ^  v2 

D  =  (cos  ft  +  U)2  +  v2  (6) 

E  =  ((m^2  _  m.2)  cos  ft  -  U)2  +  (2m^mjCos  ft  V)2 

F  =  ((m^2  .  m.2)  +  yj2  +  (2m^mjCos  ft  +  V)2 

R,2  =  C  D 
R2^  =  E  F 

The  reflectivity  is  given  by  equation  (7). 

R(ft.k)  =  0.5  (R,2  +  R22)  (7) 

For  each  wave  number  k,  the  complex  index  of  refraction  was  determined.  The  values 
used  were  taken  from  Querry  et  al.  (1977)  and  (Hale  and  Querry  (197.1)).  Figure  2  shows  the 
real  and  imaginary  part  of  the  refractivity. 
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Figure  2.  Real  and  imaginary  parts  m^.,  m  of  refractivity  as  a  function 
of  wavenumber  k. 

A  study  of  the  dependence  of  reflectivity  and  emissivity  on  frequency  showed,  that  in 
the  8-  to  I2-/Liin  region  the  variability  of  these  parameters  is  less  than  5  percent.  In  terms  of 
radiation  temperature,  the  changes  in  emissivity  and  reflectivity  could  result  in  deviations  of 
radiation  temperature  of  about  one  degree.  Since  the  effects  of  emissivity  and  reflectivity  are 
counteracting,  it  is  not  clear  what  the  overall  effect  would  be.  A  closer  look  at  this  problem  is 
necessary  because  it  lari;ely  influences  the  computation  time  of  the  model. 

2.6  NUMERICAL  EVALUATION  OF  THE  EQUATIONS 

Equation  (2)  cannot  be  solved  analytically  and  the  integration  must  be  approximated 
numerically.  The  reflection  angles  and  zenith  angles  r-ust  be  calculated  with  equations  (3)  and 
(4)  for  a  given  sensor  zenitli  angle.  The  wave  slopes  are  taken  out  of  the  intervals  30^  <Sx< 
30^.  -3ay  <  Sy  <  Oy  with  an  increment  DS.  The  sky  radiance  for  the  different  zenith  angles  is 
calculated  by  LOWiRAN  6.  To  limit  the  number  of  calculations,  the  zenith-angles  are  divided 
into  clashes.  The  results  presented  here  were  achieved  with  a  maximum  value  of  30  different 
zenith-angle  classes.  The  number  of  classes  was  further  reduced  by  the  criterion  that  each  class 
should  contain  at  least  10  percent  of  the  probability.  If  a  class  does  not  exceed  this  value,  it  will 
be  added  to  the  next  class  and  exceed  the  limit.  For  each  zenith-angle  class  i,  the  average  zenith 
angle  /a(i),  the  probability  of  this  class  P(i)  and  the  average  reflection  angle  fl(i)  is  calculated  in 
equations  8,  9,  and  10. 


P(i)  = 

SPfM)  . 

(8) 

A* 

Aifi)  = 

S/iP(M)/P(i)  .  and 

(9) 

n(i)  = 

5:n(M)P(Ai)/P(i)  . 

(10) 

Ai 
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For  the  average  zenith  angle,  the  sky  radiance  was  calculated  by  LOWTRAN  6.  This 
radiance  was  multiplied  with  the  reflectivity  for  this  zenith  angle  class  R(i)  and  with  the  proba¬ 
bility  P(i).  Subsequently,  the  values  of  all  zenith  angle  classes  were  added.  The  emitted  radiance 
of  the  ocean  was  accounted  for  by  an  averaged  frequency-dependent  emissivity  EMIS.  EMIS  is 
calculated  as  an  weighted  average  over  all  zenith  angle  classes  (equation  (11)). 

EMIS  =  SP(i)i(l  -  R(fKi)  ,  k))/  5;P(i)i  .  (11) 

The  radiance  emitted  by  the  ocean  and  the  sum  of  the  reflected  sky  radiances  were 
multiplied  by  the  transmissivity  between  sensor  and  ocean  surface  to  account  for  attenuation 
along  the  path,  then  added  to  the  emission  of  the  air  between  the  sensor  and  ocean  surface  to 
obtain  the  total  radiance. 

To  accomplish  all  the  different  functions,  several  subroutines  had  to  be  added  to  the 
LOWTRAN  6  code. 


3.0  PROGRAM  CHANGES  IN  LOWTRAN  6 
3,1  PROGRAM  STRUCTURE 

Figure  3  shows  the  structure  ot  the  LOWTRAN  6  main  program  together  with  differ¬ 
ent  subprograms.  The  structure  of  the  subroutine  GEO  is  shown  in  figure  4.  The  subroutines 
WAVE,  WAHR,  GEOPRO,  and  REFLEX  have  been  added  to  calculate  the  wave-slope  statis¬ 
tics  and  the  geometry.  The  subprogram  TRANS  is  modified  through  the  introduction  of  the 
subroutine  FRESNEL  (figure  5).  FRESNEL  calculates  the  reflectivity  and  emissivity  according 
to  equation  6. 


I  GEO  I 

I  SUB-  I 

I  ROUTINES  I 

I _ I 


(-igure  3.  LOWTRAN  6  program  structure. 
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Figure  4.  The  structure  of  subprogram  GEO  -  modrfred  by  WAVE,  WAHR, 
GEOPRO  and  REFLEX. 


Figure  5.  Subroutine  FRESNEL  Is  added  to  the  subprogram  TRANS. 


3.2  PROGRAM  FUNCTIONS 

As  few  as  possible  changes  were  made  within  the  original  LOWTRAN  code  to  allow 
the  modified  version  to  run  in  the  way  in  which  the  users  were  accustomed. 

When  the  subprogram  FNDHMN  detects  that  the  optical  path  intersects  the  earth,  the 
subroutine  WAVE  will  be  activated  to  calculate  the  wave-slope  statistics.  WAVE  creates  a 
LOWTRAN-compatible  data  set  (file:  WAVE.DAT)  of  zenith  angles  tha*  allows  radiance  to 
be  reflected  into  the  line-of-sight.  Thereafter,  the  program  returns  into  the  LOWTRAN  main 
program  to  restart  the  calculation.  The  new  data  set  will  be  read  and  the  radiances  will  be 
calculated.  In  the  subprogram  TRANS,  the  radiances  will  be  multiplied  with  the  probability  of 
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the  corresponding  zenith  angle  classes.  The  subroutine  FRESNEL  calculates  'he  icflectivity  for 
each  wave  number.  After  calculating  the  radiances  for  the  introduced  zenith  angles,  the  result 
of  equation  (2)  is  calculated  and  the  program  continues  at  the  point  where  it  was  interrupted. 

To  assure  that  all  the  different  subroutines  work  together,  several  changes  had  to 
be  made  within  the  LOWTRAN  main  program  and  in  the  subroutines  GEO,  GEOINP, 
FNDHMN,  and  TRANS.  Those  changes  were  marked  in  the  source  code  (CUUtt.  .##).  Two 
new  common  blocks,  together  with  new  control  variables  (ISEM,  IDISEM,  1ST)  were  intro¬ 
duced.  ISEM  controls  the  process  when  the  optical  path  intersects  the  earth.  1ST  controls  the 
integration  of  equation  (2)  within  the  subroutine  TRANS. 

The  variables  ISEM  and  IDISEM  were  introduced  into  the  subprograms  GEO  and 
GEOINP.  According  to  their  value,  the  program  continues  the  normal  execution  (ISEM  =  0  or 
I  respectively  ISEM  =  2  and  IDISEM  =  0)  or  it  returns  into  the  main  program  (ISEM  =  2  and 
IDISEM  >  0)  when  the  optical  path  hits  the  earth.  These  parameters  assure  a  coordinated 
program  variables. 

The  subprogram  FNDHMN  detects  whether  the  earth  is  hit  by  an  optical  path  or  not 
and  accordingly  the  value  of  the  variable  ISEM  is  set.  For  ISEM  =  I  (hit  surface),  the  subrou¬ 
tine  WAVE  is  called,  the  wave-slope  statistic  is  calculated,  and  ISEM  is  set  to  2  by  WAVE. 

The  source  code  of  the  modified  subroutines  of  LOWTRAN  6,  together  with  the  added 
subprograms  can  be  found  in  Appendix  2. 


4.0  EXPERIMENTAL  INVESTIGATION 

To  verify  the  model,  several  measurements  were  taken  to  provide  sea  and  sky  radiances 
data  together  with  meteorological  information. 

4.1  AGA  MEASUREMENTS  OF  8-  TO  I2-Mm  RADIANCE 

7  he  8-  to  l2-/im  radiance  was  measured  with  a  calibrated  thermal-imaging  system 
(AGA  Thermovision  Model  780).  The  lens  had  a  field-of-view  of  2.95°,  and  an  instantaneous 
field-of-view  of  0.9  milliradian  (mr).  The  system  response  is  determined  by  placing  a  blackbody 
of  known  temperature  (±  0.l°C  accuracy)  in  front  of  the  system  so  that  it  filled  the  instantane¬ 
ous  field-of-view.  The  digitized  video-signal  transfer  signal  of  the  system  then  allows  the  tem¬ 
perature  of  the  blackbody  to  be  reproduced  within  ±  0.2°C.  To  measure  sea  and  sky  radiance 
the  system  was  positioned  pointing  west  towards  the  ocean  on  the  Point  Loma  peninsula  in 
.San  Diego.  CA.  7  he  height  was  about  .^,7  m  above  mean  sea  level. 

4.2  MEASUREMENT  OF  METEOROLOGICAL  PARAMETERS 

The  LOWTRAN  6  model  needs  vertical  profiles  of  temperature,  relative  humidity, 
and  pressure  as  an  input  to  calculate  the  horizon  and  sky  radiances.  To  measure  the  tempera¬ 
ture,  humidity,  and  pressure,  two  different  systems  were  used.  A  Vaisala  RS  80  radiosonde, 
used  as  a  dropsonde,  was  released  from  a  Piper  Navajo  aircraft.  The  other  system  was  the 
aircraft  equipped  with  the  Rosemount  temperature  and  pressure  probes,  an  EG&G  dewpoint 
sensor,  and  a  Barnes  PRT  5  radiation  thermometer  to  measure  the  sea-surface  temperature. 
After  releasing  the  dropsonde  at  the  9000  foot  altitude,  the  aircraft  spiraled  down  to  an  alti¬ 
tude  of  100  feet.  Both  systems  gathered  a  temperature,  humidity,  and  pressure  profile  during 
the  descent. 
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Comparisons  between  the  systems  brought  up  questions  about  the  accuracy  of  the 
radiosonde  system.  The  accuracy  of  the  pressure  sensor  is  given  with  ±  0.5  hPa.  The  pressure 
height  curves  suggest  that  the  accuracy  might  be  closer  to  ±  1 .0  hPa.  The  visual  humidity  values 
were  especially  low  in  clouds.  As  a  consequence,  the  aircraft  system  was  used  in  most  cases. 

43  AVAILABLE  DATA  SETS 

The  inaccuracy  of  local  forecasts  of  fog/ stratus  dissolution  times  made  it  impossible  to 
gather  radiance  measurements,  together  with  aircraft/ radiosonde  data  for  more  than  three 
days.  Table  1  lists  the  measurement  days  and  the  systems  used. 

Table  1.  Measurement  days  and  systems. 

30  June  1987  AGA  +  Vaisala  RS  80  +  Piper  Navajo 

8  September  1987  AGA  +  Piper  Navajo 

17  September  1987  AGA  +  Piper  Navajo 

Figures  6.  7,  and  8  show  measured  temperature  and  humidity  profiles  as  a  function  of 

altitude. 
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Figure  6.  Soundings  of  temperature  and  relative  humidity  for  30  June  1987. 
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Figure  7.  Soundings  of  temperature  and  relative  humidity  for  8  September  1987. 
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Figure  8.  Soundings  of  temperature  and  relative  humidity  for  17  September  1987. 
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5.0  SIMULATION  RESULTS 


This  paragraph  presents  simulation  results  of  the  model  and,  whenever  possible, 
compares  them  with  measurements  of  the  sea  and  sky  radiances.  The  zenith  angles  are  refer¬ 
enced  to  a  sensor  height  of  33  m  above  mean  sea  level.  The  measured  temperature  and  humid¬ 
ity  profiles  were  extrapolated  up  to  40  km  using  a  standard  atmosphere. 

5.1  DEPENDENCY  OF  RADIANCE  ON  WIND  VELOCITY 

Figure  9  shows  the  radiance  calculated  with  the  meteorological  data  of  30  June  1987. 
The  Maritime  aerosol  modeT  was  used  with  the  visibility  set  to  72  km.  The  radiance  of  sea  and 
sky  is  shown  as  a  function  of  zenith  angle  for  5  different  wind  velocities  (WS  =  0.1,  1.0,  5.0, 
10.0  and  20.0  m/s). 


88  90  92  94  96  98  100 

ZENITH  ANGLE,  degrees 


Figure  9.  Radiance  as  a  function  of  zenith  angle  and  wind  velocity. 

A  common  feature  of  all  curves  is  the  steep  decrease  of  radiance  for  zenith  angles 
below  the  horizon.  As  the  zenith  angle  is  increased,  a  minimum  radiance  is  reached.  After  that, 
the  radiance  rises  again.  This  increase  in  radiance  can  be  explained  by  the  greater  importance 
of  emitted  sea  radiance.  The  behaviour  of  the  curves  in  this  region  is  determined  by  the  water 
temperature  and  the  windspeed  that  controls  the  reflectivity/ emissivity.  For  velocities  less  than 
1.0  m/s,  the  radiances  just  below  the  horizon  are  relatively  high.  The  smoother  sea  surface 
reflects  sky  radiance  just  above  the  horizon  with  relatively  high  radiation  temperature.  For 
higher  windspeeds,  sky  radiance  from  lower  zenith  angles  (lower  radiation  temperature)  is 
reflected  into  the  line-of-sight.  In  cases  of  low  windspeeds,  the  minimum  in  radiance  is  found 
at  large  zenith  angles  (98°). 
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5.2  INFLUENCE  OF  WIND  AND  VIEWING  DIRECTION  ON  RADIANCE 

The  calculated  radiance  is  plotted  for  angle  differences  between  wind  direction  and 
sensor  azimuth  as  a  function  of  zenith  angle  (figure  10).  Again,  the  meteorological  profiles  of 
30  June  1987  were  used  together  with  the  maritime  aerosol  model  and  a  visibility  of  72  km. 
The  calculations  were  for  two  different  windspeeds,  1.0  and  10.0  m/s. 


Figure  10a.  The  influence  of  wind-  and  viewing  direction  on  radiance.  WS  =  1 .0  m/s. 


Figure  10b.  The  influence  of  wind-  and  viewing  direction  on  radiance.  WS  =  10.0  m/s. 


The  radiance  was  not  affected  up  to  a  zenith  angle  of  about  91°.  As  the  zenith  angle 
increases  further,  the  radiances  differ  slightly  for  the  various  angle  differences.  The  influence  of 
this  parameter  is  not  as  important  as  the  influence  of  wind  velocity  on  radiance.  The  wind 
direction  dependency  may  be  omitted  and  the  two  dimensional  gaussian  distribution  replaced 
by  a  one  dimensional  gaussian  distribution. 

S3  CASE  STUDIES 

The  radiances  calculated  using  the  meteorological  data  of  30  June  1987,  8  September 
1987,  and  17  September  1987  will  be  compared  with  the  radiances  measured  with  the  AGA 
system.  In  all  cases,  the  Navy  Maritime  aerosol  model  (Gathman,  1983)  was  used.  The  aerosol 
size  distribution  of  this  model  is  controlled  by  the  current  windspeed,  the  24-hour  average  of 
the  windspeed  and  a  constant  airmass  (AM)  factor,  representing  the  origin  of  the  airmass.  The 
default  value  of  AM  in  LOWTRAN  6  was  used  for  all  calculations.  Wind  velocities  were 
measured  at  the  location  where  the  AGA  system  was  positioned.  Water  temperatures  were 
determined  with  the  airborne  PRT  5  radiation  thermometer.  Windspeeds  and  water  tempera¬ 
tures  are  given  in  the  legend  to  each  figure. 

Figure  1 1  shows  the  measured  as  well  as  the  calculated  radiance  for  30  June  1987.  The 
agreement  between  calculated  and  measured  radiances  was  good.  The  deviations  between  both 
curves  can  be  explained  by  the  kelp  fields  that  affected  the  measured  radiances.  The  effects  of 
the  kelp  on  sea  surface  roughness  can  be  visually  seen.  In  the  area  of  the  kelp,  the  sea  surface 
was  smoother.  Short  and  steep  waves  appeared  to  be  damped.  Aircraft  measurements  showed 
that  the  sea  surface  temperature  was  also  affected.  In  the  area  of  the  kelp,  and  between  kelp 
and  shore  the  temperature  was  higher  than  in  the  open  ocean.  But  even  with  those  deviations, 
the  agreement  between  measured  and  calculated  values  was  surprisingly  good. 


Figure  1 1.  Comparison  between  measured  and  calculated  radiances  for 
30  June  1987. 


The  data  of  8  September  1987  can  be  compared  with  30  June  1987  in  respect  to  winds- 
peed  and  water  temperature.  But,  as  figure  12  shows,  the  agreement  between  calculated  and 
measured  radiances  was  better  on  8  September  1987.  On  this  day,  the  visibility  was  very  poor 
—  the  Coronado  islands,  about  30  km  away,  were  not  visible.  This  might  have  contributed  to 
the  better  agreement  since  we  expected  a  stronger  influence  of  the  atmosphere  between  sensor 
and  sea  surface  and  less  contribution  by  the  ocean. 


Figure  12.  Measures  and  calculated  radiances  as  a  function  of  zenith  angle 
for  8  September  1987. 


The  results  of  17  September  1987  are  presented  in  figure  13.  This  data  show  the  largest 
differences  between  calculated  and  measured  radiances  (the  differences  in  radiance  are  propor¬ 
tional  to  differences  of  less  than  1.0°C  in  radiation  temperature).  There  are  several  possible 
explanations.  First  of  all,  the  temperature  and  humidity  profiles  were  available  only  up  to  a 
height  of  1.3  km  whereas  in  the  other  cases  the  profiles  reached  up  to  2.8  km.  That  affected  the 
calculations  (we  do  not  know  how  much).  A  second  possibility  is  that  the  influence  of  kelp 
fields  changed  with  windspeed. 
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Figure  13.  Measured  and  calculated  radiances  as  a  function  of  zenith  angle 
for  17  September  1987. 


6.0  CONCLUSIONS 

The  paper  showed  that  the  suggested  sea-radiance  model  was  able  to  predict  sea 
radiances  both  qualitatively  and  quantitatively  correct  in  three  cases.  Further  studies  should 
show  how  often  a  simple  sea-state  model  can  be  used  and  under  which  circumstances  it  fails. 

Further  investigations  should  address  the  maximum  height  requirement  at  which  the 
temperature  and  humidity  profiles  need  to  be  measured  to  obtain  better  agreement  between 
measured  and  calculated  radiances. 

The  maximum  number  of  required  zenith  angles  representing  the  measured  radiances 
needs  to  be  determined.  Also  the  validity  of  a  single  frequency-dependent  emissivity  and  reflec¬ 
tivity  needs  to  be  investigated. 
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7.0  SYMBOLS  AND  ABBREVIATIONS 


BB(k,Ta) 

EMIS 

k 

P(m) 

P(Sx,Sy) 

P(i) 


Po 

R(ak) 

Rl 

R2 

S(0,k) 

s(e,k) 

Sv 


S| 

Soo 

Ta 

Tw 

a 

p 

0 

0’ 


/x(i) 

o'x 

ay 

T 

’’a 

’"s 

n 

fKi) 

IR 

P 

RH 

WS 


Planck  blackbody  function 
mean  emissivity  of  the  sea 
wave  number 

probability  that  radiance  from  the  zenith  angle  n  reaches 
the  sensor 

probability  of  the  wave  tilts  Sy 

probability  that  radiance  from  the  i-th  zenith  angle  class 

reaches  the  sensor 

probability  for  shadowing 

reflectivity 

reflectivity  for  parallel  polarized  waves 
reflectivity  for  perpendicular  polarized  waves 
radiance 
mean  radiance 

wave  tilt  in  x-direction  =  tan  a 
wave  tilt  in  y-direction  =  tan 
position  on  the  optical  path 
sensor  position 
reflector  position 
infinity 

air  temperature 

water  temperature 

angle  of  wave  tilt  in  x-direction 

angle  of  wave  tilt  in  y-direction 

sensor  zenith  angle 

sensor  zenith  angle  at  the  reflection  point  modified  by 
atmospheric  refraction 

zenith  angle  from  where  radiance  is  reflected  into  the 

observer  line  of  sight 

mean  zenith  angle  of  the  i-th  class 

standard  deviation  of  the  wave  tilt 

standard  deviation  of  the  wave  tilt  Sy 

total  transmission  = 

absorption  transmission 

scattering  transmission 

angle  difference  between  wind  direction  and  sensor 

azimuth 

reflection  angle 

mean  reflection  angle  of  the  i-th  class 

infrared 

pressure 

relative  humidity 
windspeed 
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APPENDIX  A:  SHADOWING  OF  WAVES 


The  shadowing  of  waves  by  other  waves  is  a  problem,  especially  for  zenith  angles  close 
to  the  horizon.  Other  authors  limit  the  possible  zenith  angle  range  of  their  models  to  values  far 
away  from  the  horizon  (Sidran,  1981;  Schwartz  and  Hon,  1986)  or  they  omit  the  contributions 
entirely  (Saunders,  1968). 

In  this  paper,  the  waves  are  assumed  to  be  statistically  independent  from  one  another. 
Under  this  assumption,  McAdam  (1969)  (in  the  paper  of  Gordon,  1969)  showed  that  the  prob¬ 
ability  for  shadowing  depends  only  on  the  incident  angle  (relative  to  the  vertical)  and  on  the 
standard  deviation  of  the  wave  slopes.  If  we  assume  that  the  wave  slope  distribution  for  the 
shadowing  waves  is  also  given  by  the  Cox-Munk  wave-slope  statistic,  the  shadowing  effects  can 
be  determined  by  the  probabilities  for  the  occurrence  of  the  zenith  angle  /i. 

For  a  given  set  of  parameters,  0,  ax,  ay,  the  probability  for  shadowing  is  Pq.  In  this 
case,  the  optical  paths  contained  in  Pq  do  not  hit  the  wave  facet  under  consideration.  The 
optical  paths  hit  a  facet  out  of  a  normally  distributed  ensemble.  These  rays  are  again  shadowed 
with  Pq  and  so  forth.  When  P(/z)  is  the  probability  (including  shadowing)  that  radiance  is 
reflected  from  the  zenith  angle  n  into  the  line-of-sight,  we  have  the  probability  P|(/i)  =  (1  -  Pq) 
P(^).  Of  the  shadowed  optical  paths,  a  percentage  P2  =  Po(I  -  Po)P(m)  of  radiance  is  reflected 
from  the  zenith  angle  m  into  the  sensor.  This  can  be  repeated. 

1.  shadowing  Pj(m)  =  (1  Po)*P(m) 

2.  shadowing  P2(m)  =  Po*(l  Po)*P(m) 

3.  shadowing  PjfM)  =  Po^*(l  -  Po)*P(m) 


n.  shadowing  Pn<M)  =  Po"*(l  Po)*P(m) 

If  the  single  contributions  are  added,  we  get  for  n  —  00; 

Sn  =  (l  Po")P(M).  Po<l 

lim  I  -  Pq"  =  I 

Sn  =  P(M) 

After  a  large  number  of  steps,  we  obtain  the  true  probability  P(/u)  for  radiance  reflected 
from  this  zenith  angle.  It  is  possible  to  simplify  the  calculation  by  integrating  over  the  zenith 
angle  after  the  first  step  to  obtain 

P(^)  =  p,(^)/(i  -Pq)  =  p,/5:p,  . 
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